Fluctuoscopy of Disordered Two-Dimensional Superconductors 



A. Glatz, 1 A. A. Varlamov, 1 ' 2 and V. M. Vinokur 1 

1 Materials Science Division, Argonne National Laboratory, 

9700 S.Cass Avenue, Argonne, Illinois 60637, USA 
2 CNR-SPIN, Viale del Politecnico 1, 1-00133 Rome, Italy 
(Dated: December 7, 2010) 

We revise the long studied problem of fluctuation conductivity (FC) in disordered two-dimensional 
superconductors placed in a perpendicular magnetic field by finally deriving the complete solution 
in the temperature-magnetic field phase diagram. The obtained expressions allow both to perform 
straightforward (numerical) calculation of the FC surface ScTxx (T,H) and to get asymptotic ex- 
pressions in all its qualitatively different domains. This surface becomes in particular non-trivial 
at low temperatures, where it is trough-shaped with 5oxx (T,H) < 0. In this region, close to the 
quantum phase transition, Sa^x^iT, H — const) is non-monotonic, in agreement with experimental 
findings. We reanalyzed and present comparisons to several experimental measurements. Based 
on our results we derive a qualitative picture of superconducting fluctuations close to H C 2 (0) and 
T = where fluctuation Cooper pairs rotate with cyclotron frequency iv c ~ AgJ s and Larmor 
radius ~ £bcs , forming some kind of quantum liquid with long coherence length £qf 3> £bcs and 
slow relaxation (tqf 3> fiAgJ s ). 



I. INTRODUCTION 

The understanding of the mechanisms of supercon- 
ducting!-, fluctuations (SF), achieved during the past 
decadcsu provided a unique tool obtaining information 
about the microscopic parameters of superconductors 
(SC). SFs are comprised of Cooper pairs with finite life- 
time which appear already above the transition but do 
not form a stable condensate yet. They affect thermody- 
namic and transport properties of the normal state both 
directly and through the changea.wh.ich they cause in the 
normal quasi-particlc subsystem^. 

SFs are commonly described in terms of three princi- 
pal contributions: the Aslamazov-Larkin (AL) process, 
corresponding tp the opening of a new channel for the 
charge transfer! 2 ], anomalous Maki- Thompson (MT) pro- 
cess, which describes single-particle quantum interference 
at impurities in the presence of SFsoQ, and the change 
of the single-particle density of states. (DOS) due to their 
involvement in fluctuation pairingstM. The first two pro- 
cesses (AL and MT) result in the appearance of positive 
and singular contributions to conductivity (diagrams 1 
and 2 in Fig. Q) close to the superconducting critical 
temperature T c0 , while the third one (DOS) results in 
a decrease of the Drude conductivity due to the lack of 
single-particle excitations at the Fermi level (diagrams 
3-6 in Fig. ^). The latter contribution is less singular 
in temperature than the first two and can compete with 
them only if the AL and MT processes are suppressed 
for some reasons (for example, c-axis transport in lay- 
ered superconductors) or far away from T c q. 

The classical results obtained first in the vicinity of T c q 
were later generalized to temperatures far from the tran- 
sition, e.g. in Refs. and to relatively high fields, 
see Ref. |Uj. More recently, quantum fluctuations (QF) 
came into the focus of investigations. In Ref. |ll] and Q 
it was found, that in granular SC at very low tempera- 
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FIG. 1. (Color online) Feynman diagrams for the leading- 
order contributions to the electromagnetic response operator. 
Wavy lines stand for fluctuation propagators, solid lines with 
arrows are impurity-averaged normal state Green's functions, 
crossed circles are electric field vertices, dashed lines with a 
circle represent additional impurity renormalizations, and tri- 
angles and dotted rectangles are impurity ladders accounting 
for the electron scattering at impurities (Cooperons). 



tures and close to H C 2 (0), the positive AL contribution to 
magneto-conductivity decays as T 2 while the fluctuation 
suppression of the DOS results in a temperature inde- 
pendent negative contribution, logaritmically growing in 
magnitude for H — > H C 2 (0). The authors of Ref. IDj] 
came to the same conclusion while studying the effect 
of QFs on the Nernst-Ettingshausen coefficient in two- 
dimensional (£p) SCs. For the first time they attracted 
the attentionliil to the special role of diagrams 9 & 10 in 
Fig. | 

Special attention should be paid to the paper by Gal- 
itski and Larkin, Ref. |l5|, where the effects of QFs 
on magneto-conductivity and magnetization of 2D su- 
perconductors were studied. The authors analyzed all 
ten diagrams shown in Fig. [l] in the lowest Landau level 
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FIG. 2. (Color online) Fluctuation correction to conductivity 



(FC) Sa 



(tot) 



(t, h) as function of the reduced tempera- 



ture t = T/T c0 and magnetic field h = 0.69H/H c2 (0) plot- 
ted as surface. The FC changes its sign along the thick red 
line (Sa = 0). The boundary of the superconducting region 
is shown by a dashed line. Here Sa is plotted for constant 
tT c0 = 1CT 2 and t^Tcq = 10. 



(LLL) approximation, valid at fields close to the criti- 
cal line H C 2 (T). They found an expression for the total 
fluctuation contribution to magneto-conductivity in this 
region and demonstrated that, analogously to the situ- 
ation in granular SCs, close to zero temperature and in 
the vicinity of H C 2 (0) this contribution is negative, i.e. 
QFs increase resistivity, and not conductivity (in con- 
trast to the situation close to T c p). Nevertheless, con- 
trary to the conclusions of Refs. Jll] and [l^] , in Rcf. 
the logarithmic growth of the magneto-resistance at zero 
temperature when the magnetic field approaches H C 2 (0) 
from larger fields, is due to all AL, MT and DOS-like 
contributions. 

Yet, the confidence in the exotic nature of negative 
fluctuation corrections and the common believe of fluc- 
tuation contributions to conductivity being positive be- 
yond the narrow domain of the quantum phase transition, 
has been persistent and is based on available asymptotic 
expressions only. The region near T = and magnetic 
fields near H C 2 (0) remains poorly understood and in ad- 
dition, an universal picture combining QFs at high mag- 
netic fields and conventional finite temperature quantum 
corrections is still lacking. 

This is why we revisit the problem of fluctuation con- 
ductivity of a disordered 2D superconductor placed in a 
perpendicular magnetic field in this paper. We present 
an exact calculation (without the use of the LLL approx- 
imation) of all ten diagrams of the first order of fluctu- 
ation theory (see Fig. [j]) valid in the whole H-T phase 
diagram beyond the superconducting region, i.e. for ar- 
bitrary fields H > H C 2 (T) and temperatures T c (H) < T. 
The obtained expressions allow both to perform straight- 
forward (numerical) calculation of the fluctuation con- 
ductivity "surface" Saxx (T, H) and to get asymptotic 
expressions in all its qualitatively different domains. 

A typical example of the surface Saix (T, H) is prc- 
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FIG. 3. (Color online) Contours of constant fluctuation con- 
ductivity [8a = Sax°^ (t, h) shown in units of e 2 ]. The domi- 
nant FC contributions are indicated by bold-italic labels. The 
dashed line separates the domain of quantum fluctuations 
(QF) [dark area of Sa < 0] and thermal fluctuations (TF). 
The contour lines are obtained from Eq. (^|) with T c qt = 0.01 
and T c0 "i> = 10. 



scntcd in Fig. g and demonstrates that our revision and 
completion of the commonly believed understanding of 
fluctuation corrections is urgently called for: Its strik- 
ing feature consists of the fact, that the FC is positive 
only in the domain bound by the separatices H C 2(T) 

and Sa^i^H) = and is negative throughout all 
other parts of the phase diagram (see Fig. ||, in which 

the domains of different overall signs of 8<Jxx (T, H) 

and contours of constant Saxx in the whole phase di- 
agram are shown). Contrary to the common assump- 
tion, the FC is only positive in the domain of weak 
fields and temperatures above T c o, the region of posi- 
tive corrections depends on the magnitude of the posi- 
tive anomalous MT contribution (i.e. on the value of the 
phase-breaking time r^). With increasing magnetic field, 

the interval of temperatures where 5axx^ (T, H) > 
shrinks and becomes zero close to H C 2(0). In particu- 
lar at low temperatures, the behavior of the FC turns 
out to be highly non-trivial. In this case, the surface 



S&xx (T, H) has a trough-shaped character and the de- 
pendence 5aix(T,H = const) is non-monotonic. We 
will see below that this feature is observed in available 
experimental results as well. 

Our analysis also elucidates the understanding of the 
hierarchy of the various contributions to the fluctuation 
corrections in different domains of the phase diagram (see 
Fig. U in which the dominating fluctuation contributions 
to magneto-conductivity are indicated for different re- 
gions of the phase diagram). We reproduce the expres- 
sion for the total fluctuation contribution to magncto- 
conducitivty of 2D disordered superconductors presented 
in Rcf. [15 in the vicinity of the H C 2(T) line. Neverthe- 
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FIG. 4. (Color online) Illustration of the cluster structure of 
a FCP (2e) liquid above the upper critical field. This picture 
represents a snapshot at a certain time and would stay that 
way for time tqf- The typical size of a coherent FCP cluster 
is £qf- 



less, our analysis clearly demonstrates that the princi- 
pal fluctuation contributions close to T c q, paraconduc 
tivity (AL), anomalous MT and DOS, in the region of 
QF become zero as ~ T 2 (compare to Refs. jy] and 
- which is contrary to the picture given in Rcf. Jla 
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It is 

the fourth, usually ignored, fluctuation contribution, for- 
mally determined by the sum of diagrams 7-10 and the 
regular part of the MT diagram, which governs the quan- 
tum phase transition (QPT). It can be identified with 
the renormalization of the single-particle diffusion coeffi- 
cient in the presence of fluctuations (DCR) and it turns 
out that this contribution dominates in the periphery 
of the phase diagram including the vicinity of the QPT 

t = T/T c0 <h=(H- H c2 (0)) /H c2 {0):H > H c2 (0) . 

Finally, based on our results, we propose a qualita- 
tive picture for QPT, which drastically differs from the 
Ginzburg-Landau one, valid close to T c o- The latter 
can be described in terms of a set of long-wavelength 
fluctuation modes [with A > £gl (T) S> £bcs] of the 
order parameter, with characteristic lifetime tql = 
irh/8k B (T - T c0 ). Near the QPT, the order parame- 
ter oscillates on much smaller scales - the fluctuation 
modes with wave-lengths up £bcs are excited. Due to 
the magnetic field, one can imagine that FCP in this re- 
gion rotate with the Larmor radius ~ £bcs and cyclotron 
frequency lo c ~ ^bcs- We snow that close to H c2 (0) 
these FCP form some kind of quantum liquid with long 
coherence length £qf ~ ^Bcs/h 1 ^ 2 and slow relaxation 
tqf ~ hA^/h (see Fig. §. 

In the following sections and in the appendices, we 
will show the details of our derivations and calculation 
and present the general expression for the fluctuation 
magneto-conductivity of disordered 2D SC throughout 
the whole phase diagram. For the calculation of the com- 
plete and various fluctuation corrections (by numerical 
integration and summation), we developed an optimized 
program which is available at |l6| . It can be used as the 
theoretical basis for the "fluctuoscopy" of superconduc- 
tors: the study of their behavior in ultra-high magnetic 
fields and precise extraction for their physical parame- 
ters like the critical temperature and magnetic field, and 



the temperature dependence of the phase-breaking time, 
and/or e.g. for the separation of the quantum corrections 
in studies of the "superconductor-insulator" transition. 



II. MODEL 

We consider a disordered 2D superconductor charac- 
terized by the diffusion coefficient T> placed in a perpen- 
dicular magnetic field H at temperatures T > T c (H) . 
Temperatures should not be too close to the critical tem- 
perature and remain beyond the region of critical fluctu- 
ations, i.e. T/T c (H) - 1 > y/Gi^jH). The Ginzburg- 
Levanyuk number Gi(2) for conductivity (sec Rcf. 
in both extremes of the line H c2 (T)(aX temperatures 
close to T c o and at zero temperature) is on the order 
of {p 2 F ld) , where d is the SC film thickness, and it can 
reach values of up to 10 ~ 2 . We assume the temperature 
T <C min {t -1 ,uid} in order to remain in the diffusive 
regime of electron scattering and in the frameworks of 
the BCS model (r is the electron elastic scattering time 
on impurities, wd is the Dcbcy frequency). The restric- 
tions on magnetic field are dictated by the requirements 
to be below the regime of Shubnikov-de Haas oscillations 
[w c t < 1 H < (r c0 r) _1 H c2 (0) , where w c = WeH 
is the fluctuation Cooper pair cyclotron frequency] and 
to be below the Clogston limit: H < (eft) H c2 (0) , i.e. 

H/H c2 (0) « min {(Tcor)" 1 ,e F r}. 

Under these rather non-restrictive assumptions the DC 
fluctuation conductivity 



<5cr (fl) (T, H) = - lim 



Ini 



(1) 



is determined by the imaginary part of the fluctuation 
contribution Qw(w,T, H) to the electromagnetic re- 
sponse opcratoiu. The latter is described graphically by 
the ten standard diagrams shown in Fig. |l|. The solid 
lines denote the one-electron Green function 
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wavy lines correspond to the fluctuation propagator 
L- 1 (V k )= (2) 

\n k \ +uj c (n + iy 
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and shaded three- and four-leg blocks indicate the results 
of the average over elastic impurity scattering of electrons 
(Cooperons): 



A„ (si,e 2 ) 



r-^(- ei e 2 ) 
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Here i/q is the one-electron density of states, n, m are 
the quantum numbers of the Cooper pair Landau states, 
Q k = 2nkT , ei = 2nT(l + l/2) arc the bosonic and 
fermionic Matsubara frequencies. An important charac- 
teristic of these expressions is that they are valid even 
far from the critical temperature [for temperatures T <C 
min{T~ , W£>}] and for <C w.d an d n -C (T c ot) 

In the Appendices we present the details of the calcula- 
tion of all ten diagrams performed under the above gen- 
eral assumptions. In the following sections of the main 
text, we restrict ourselves to the discussion and analysis 



of the main result: the complete expression of the fluctu- 
ations corrections and the individual contributions from 
AL, MT, DOS, and DCR processes. 



III. RESULTS 



The complete expression for the total fluctuation cor- 
rection to conductivity ScTxx (T, H) of a disordered 2D 
SC in a perpendicular magnetic field that holds in the 
complete T-H phase dia gram above the line H r o(T) is 
given by the sum of Eqs. ( Ajjh, ©, @, and (Dig): 
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Here t = T/T c0 , 



h = 



H 



8je H c2 (0) 



= 0.69 



H 



H c2 (0) : 



7_e = e 7e (7e is the Eulcr constant), M = (tT C QT) 1 , 7^ 
7r/ (8T c or^), is the phase-breaking time, 



£ m = £ m (t, h, z) = \a.t+%jj 



and its derivatives £m {t,h,z) = d^£ m (t,h,z). Apart 
from the detailed derivation of the result, Eq. (||), one 
can also do a careful study of the asymptotic expressions 
for different fluctuation contributions throughout the h- 
t phase diagram, presented in the Appendices. All of 
them, side by side with the asymptotic expressions for 
Scr^x are summarized in table |. 

We start the discussion of table | for domains I -III, 
corresponding to the Ginzburg-Landau region of fluctu- 
ations close to T c q and in zero magnetic field (domain 
I). One can see, that our general expression Eq. (||) 
naturally reproduces the well known AL, MT and DOS 
contributions. The only new result here is the explicitly 
written contribution 5a^ DCR ^ (diagrams 7-10), which was 
usually ignored in view of the lack of its divergence close 
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FIG. 5. (Color online) Schematic representation of the regions 
of different behavior of fluctuation conductivity in the h-t 
diagram. See table [o] for more explanations on the domains. 



to T c o. Nevertheless, one can see that its constant con- 
tribution ~ lnln(T c oT) 1 is necessary for matching the 
GL results with the neighboring domains VIII & IX. The 
domains II & III are still described by the GL theory in 
weak magnetic fields and Eq. (||) reproduces all available 
asymptotic expressions found in literature. 
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TABLE I. Asymptotic expressions in different domains, shown in Fig. |5[ The first column gives the domain according to that 
figure and is determined by the t & h regions given in table |n| 



domain 


t and h range 


I 


h = 0, 6 < 1 


II 


h — h c 2 ~ e <C 1 


III 


ft - 7lc2(t) < 1, e < 1 


I - III 


h < 1, e < 1 


IV 


t < ft- ftc2(t) 


V 


ft ~ t < 1 


VI 


ft < * < 1 


VII 


h - h c2 {t) < t <s: h c2 (t) 


VIII 


hit > l,h< t 


IX 


h 2> max{l, f} 



description 



zero field, near T c o 
near T^o-reflected ft C 2-line 
near /i C 2-line 
full GL region 
region of QFs 
quantum-to-classical 
classical, near h C 2(t — > 0) 
classical, strong fields 
high temperatures 
high magnetic fields 



TABLE II. Explanation of the different domains with t & h 
ranges. Here e = Int. 



The most surprising result in tabic [| is the domain IV, 
the region of quantum fluctuations (see Fig. |^): one sees 
that the positive AL (the anomalous MT contribution is 
equal to the AL one in that domain) decays with decreas- 
ing temperature as T 2 . Moreover, it is exactly cancelled 
by the negative contribution of the four DOS-like dia- 
grams 3-6: 



Sa^ = 8a?JW = -6a 



DOS 



4e 2 7 |t 2 



(4) 



tot) 

XX 



The total fluctuation contribution to conductivity 5a. 
in this important region (t <C h) is completely deter- 
mined by the renormalization of the diffusion coefficient 
(the regular part of the MT contribution and diagrams 
7-10). It turns out to be negative and at zero temper- 
ature diverges logarithmically when the magnetic field 
approaches _ff c2 (0). The non-trivial fact following from 
Eq. (§) is that an increase of temperature at a fixed value 
of the magnetic field in this domain mainly results in a 



further decrease of conductivity 

67ee 2 t 



Sa^ 



2e 2 , 1 
3vH h 



h 



O 



(5) 



and only at the boundary with domain V, when t ~ h, 
the total fluctuation contribution Saxx passes through 
a minimum and starts to grow. Such non-monotonic be- 
havior of the conductivity closc.lt) H C 2 (0) was multiple 
times observed in experimcntsLDO (see Fig. ^). 

The domain V describes the transition regime between 
quantum and classical fluctuations, while in the domains 
VI-VII, extended along the line H C 2 (T), superconducting 
fluctuations have already classical (but non-Ginzburg- 
Landau) character. In all these three regions one ob- 
serves exactly the same cancellation of the AL and DOS 
contributions as in domain IV and Saxx^ is determined 
by the negative DCR contribution. 

Finally, in the peripheral domains VIII-IX, the di- 
rect positive contribution of fluctuation Cooper pairs 
(AL) to conductivity decays faster than all the other: 
~ In -3 (T/T c o). We stress, that this exact result dif- 
fers from the evaluation of the AL paraconductivity 
far from the transition of Ref. 0, but is in complete 
agreement with the high temperature asymptotic ex- 
pression for the paraconductivity of a clean 2D super- 
conductor, see Ref. |jT9|| . This agreement seems natu- 
ral: fluctuation Cooper pair transport is insensitive to 
impurity scattering. The anomalous MT contribution, 
in complete accordance with Refs. |?] and ||, decays as 
~ ln7^V ln ~ 2 ( T / T co)- The contribution of dia grams 6- 
6 also decays as ln~ 2 (T/T c q), but without the large fac- 
tor ln7 ( ^ 1 . Finally, the regular MT contribution together 
with the ones from diagrams 7-10 decay extremely slow, 
in fact double logarithmically: 



3tt 2 I T c0 t 



T N 

— ln ln 



(G) 
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FIG. 6. (Color online) Total fluctuation conductivity for different T c qt^ for several constant temperatures t and magnetic 
fields h. a) to e) show 5a(t) for different magnetic fields below (h = 0.01,0.35, superconducting region marked by "SC"), at 
(h = 0.69), and above (h — 0.75, 1.0) the zero temperature critical field H C 2(0). The legend key of a) applies to all panels. Note, 
that the value of T c §t§ near the transition at t — can determine the sign of the fluctuation conductivity - however at very 
low temperatures the FC becomes independent of the phase-breaking time and all lines coalesce which is not resolved in that 
plot, f) to j) show Sa(h) for constant temperatures, below (t = 0.01,0.1,0.5), at (t = 1.0), and above (t — 1.5) the transition 
temperature T c o- All plots are calculated for T c ot = 10 -3 . If 5a = is within plot range, it is marked by a horizontal dashed 
line, and the critical magnetic fields H C 2(T) are shown as (red) vertical dashed lines. See detailed discussion in the text. 
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FIG. 7. (Color online) Fluctuation conductivity contribu- 
tions: AL, MT, DOS, DCR, and total (tot) for T c0 r = 10" 3 
and TcoTcj, — 5. a) shows the temperature dependence at low 
field h = 0.01, and b) the field dependence at low temperature 
t = 0.01. 



Up to the numerical prefactor this expression coincides 
with the results of Ref . and |J . 

Eq. (||) provides the basis for a "fluctuoscope" for 
superconductors, i.e. the extraction of its microscopic 
parameters from the analysis of fluctuation corrections. 
Indeed one can sec that ScTxx depends on two super- 
conducting parameters: T c q, H C 2 (0), the elastic scatter- 
ing time r, and (temperature dependent) phase-breaking 
time t<£ (T) . The elastic scattering time can be obtained 
from the normal state properties of the superconductor, 
while the Eq. (§) can become the instrument for the 
precise determination of the critical temperature T c q (in- 
stead of the often used rule "half width of transition" ) 
and H C 2 (0). Moreover, it can be an invaluable tool for 
the study of the temperature dependence of the phase- 
breaking time 70 (T). 

The exemplary surface of 8<7xx {T, H) presented in 
Fig. | for T c0 t = 10" 2 and T c0 7> = 10 shows that 
the value of determines the behavior of fluctuation 
corrections only in the region of low fields. It is con- 
venient to analyze Fig. || side-by-side with Fig. || where 
lines 8(Txx(T,H) = const through the phase diagram 
are shown. It is interesting to note that the numerical 
analysis of Eq. (^) shows that the logarithmic asymp- 



totic Eq. (^j) is valid only within an extremely narrow 
field range h < 1CT 6 . 

In order to get a broader overview of the richness of our 
main result, wc complied several magneto-conductivity 
single- parameter dependencies (cuts through the "sur- 
face" at constant t or h) in Fig. Q Each individual panel 
a) through j) of this figure shows the total FC for six 
different values of Tcqt^ between 0.1 and 50 and fixed 
T c qt = 10~ 3 . If da = is within plot range, it is 
marked by a horizontal dashed line, and the critical mag- 
netic fields H C 2(T) are shown as (red) vertical dashed 
lines, a) to e) show 6cr(t) for different magnetic fields 
below (h = 0.01,0.35, superconducting region marked 
by "SC"), at (h = 0.69), and above (h = 0.75, 1.0) the 
zero temperature critical field H C 2(0). The legend key 
of a) applies to all panels. The behavior is as expected 
from the above discussion, but it is very educative to 
take a closer look at the behavior near the QPT [panel 

c) ]: As mentioned above, the asymptotic expression for 
the quantum regime is only valid at extremely small tem- 
peratures, which cannot be resolved in this plot. There- 
fore one sees in particular for the smallest T^oT^-value 
a sharp dip in the FC at low temperatures, which will 
eventually coalesce with all other curves at even lower 
temperature (not visible) and become independent of the 
phase-breaking time, as can be seen at larger h in panel 

d) and e). 

f) to j) show Sa(h) for constant temperatures, below 
(t = 0.01,0.1,0.5), at (t = 1.0), and above (t = 1.5) 
the transition temperature T c q. Here again, it is seen 
in panel f) that in the quantum regime the FC is mostly 
independent of the phase-breaking time and only close to 
the QPT a separation of the curves becomes visible since 
even the small temperature t = 0.01 becomes of order 
h or even larger and the asymptotic expression does not 
hold anymore. 

In Fig. we plotted two particular curves of Fig. |^ in 
more detail, showing the different contributions from the 
diagram groups a)-d) of Fig. [l]. These are the curves for 
T c oT<j> = 5 at lowest magnetic field h = 0.01 [in a)] and 
temperature t = 0.01 [in bV . Comparing these curves 
to the asymptotics of table \. one sees, that the behav- 
ior near T c q is as expected [see a)] , and in particular the 
contribution from diagrams 7-10 is negligible. However, 
in the quantum regime it becomes the dominating con- 
tribution, rendering the total FC negative and only close 
to the QPT is cancelled by the MT contribution. 

Despite Eq. (||) being a closed expression, its spe- 
cific evaluation in the most general case requires sophisti- 
cated numerical summation and integration. While being 
straight-forward, one might encounter technical difficul- 
ties in the evaluation of the complex Poly- Gamma func- 
tions (z). Moreover, the summation cut-off param- 
eter Mean reach extremely large values at low tempera- 
tures [ experimental values (r c or) exp for materials near 
the superconductor-insulator transition can be on the or- 
der 10 6 ], which slows down the numerical procedure sig- 
nificantly. The latter difficulty can be partially overcome 



8 



by evaluation of the slowly divergent tails of the m-sums 
in Eq. (|J) as integrals. Here, we should also note that for 
fitting purposes one does not need to choose the real, of- 
ten extremelly small, experimental values (TcO T ) cxp - To 
save CPU time, one can assume the value (T c oT) num of 
this parameter to be much larger than (T , c o r ) cxp (but still 
much less than T^t^) and only at the very end to shift 
the final expression by In In \° a ( num . Nevertheless, the 

numerics of the problem remains challenging: for the sur- 
face plot in Fig. U we evaluated 10 6 values for 8a with 
the modest assumption (T c ot) = 0.01, yet it still took 
three month of single CPU time for its calculation. Our 
optimized tool for the evaluation of Eq. (|J) can be found 
at @. 

IV. COMPARISON WITH EXPERIMENTAL 
RESULTS 



and critical magnetic field i? C 2(0) = 13T. The resistivity 
was measured, depending on magnetic field, for low tem- 
perature values T = 200, 300, 400, 500mK. We plotted 
the theoretical expression for Sa^x^ using the fitting pa- 
rameter values B C 2(0) = 13. 7T, T c qt$ = 5 ± 1, and the 
experimentally found value of T c q = 3.35K. Overall the 
fitted FC curves show good agreement with the results 
of the measurements. 

Second we re-analyzed data for thin films of the com- 
pound La2- x Ce a; Cu04, with x = 0.09, in high magnetic 
fields, see Ref. jT^ja, with thickness lOOnm. In contrast 
to the indium oxide film, we fitted the temperature de- 
pendence of the FC at different constant high magnetic 
fields here. This sample has T c o « 22. 5K and we fitted 
B C 2(0) ~ 6T, which the 4T data suggests. In Fig. ^ we 
plotted the temperature dependence of Eq. (§) for differ- 
ent magnetic fields for fixed Tcqt^ = 0.5 and T c qt = 10~ 3 . 



A main aspect of this work is, that the complete ex- 
pression, Eq. |^, can be used to extract experimental pa- 
rameters of thin superconducting films from measured 
data ("fluctuoscopy" ). In particular the critical tem- 
perature T c0 , the critical magnetic field i/ c2 (0), and the 
phase-breaking time r^. 




exp, 500mK 

fitfor500mK 
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exp, 300mK 
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exp, 200mK 
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FIG. 8. (Color online) Comparison to resistivity measure- 
ments in thin indium oxide films, published in Ref. |^] . Here 
we present the data taken from Fig. 4a of ]2^] for the "Weak" 
sample with thickness 30nm, T c q — 3.35K, and B C 2(0) = 13T. 
We fitted the resistivity R for temperatures 0.2, 0.3, 0.4, and 
0.5K using our full expression for 5a with the experimentally 
found T c q. For B C 2(0) we fitted a slightly larger value of 13. 7T 
and T c oT^ = 5 ± 1. 

As an example of the practical use for our results, 
we fitted two different sets of experimental data. First 
we compared our general Eq. (||) to resistivity measure- 
ments in thin disordered indium oxide films, presented in 
Ref. p0[ |. Figure || shows the low temperature data for 
one sample (referred to as "Weak" in Ref. f20fl ) of a film 
with thickness 30nm, transition temperature T c q — 3.35K 
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FIG. 9. (Color online) Temperature dependence of the FC 
at different fields close to H C 2(0) and comparison to experi- 
mental data for a thin film of La2-:rCea:Cu04 (x — 0.09 ) w ith 
T c0 ~ 22. 5K and fitted£ c2 (0) w 6T (see Fig. la of Ref. Qa). 
For the theoretical curves a fixed T c ot> = 0.5 is used which 
nevertheless, captured the overall behavior. All curves are 
numerically calculated with T c qt = 10 -3 . The highest field 
curve could be possible fitted better, using a smaller r^. How- 
ever, as described in the text, is in general temperature 
dependent and the available data does not allow to extract 
this values. 



At this point it is important to remark, that de- 
pends on temperature in general, such that for a better 
fit one needs first to analyze FC data at constant temper- 
atures to extract t^(T) and then fit temperature depen- 
dent data. This way one can obtain precise values for the 
otherwise difficult to determine experimental parameters 
T c o, H c2 {0), and i>(T). However, the data of Ref. |f|a 
does not allow to extract this information for a better fit 
in Fig. | . 
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V. QUANTUM LIQUID OF FCP IN THE 
VICINITY OF H c2 (0). 

An analysis of the obtained results allows us to of- 
fer a qualitative picture of the quantum phase transition 
(QPT) occurring in the vicinity of H C 2(0) at very low 
temperatures. Above we presented the complete micro- 
scopic calculation. However, it is instructive to start our 
discussion of the QFs by describing and refreshing the 
qualitative picture of SEs in the vicinity of T c0 , in the 
Ginzburg-Landau rcgioncl, for further comparison. In do- 
mains I- III, the lifetime of fluctuation-induced Cooper 
pairs t gl can be obtained in the simplest way by using 
the uncertainty principle. Indeed, tgl ~ h/AE, where 
AE is the energy difference fcs(T — T c o) ensuring that 
tgl should become infinite at the transition point. This 
yields the standard Ginzburg-Landau time 



tgl ~ h/k B (T - T c0 ) ~ hj (k B T c0 e) 



(7) 



where e = (T — T c o) /T c q < 1 is the reduced tempera- 
ture. In its turn the coherence length £gl (T) can be es- 
timated as the distance, which two electrons move apart 
during the GL time: 

£gl (e) = (2?T GL ) 1/2 ~eBCs/V^ 



Here £bcs ~ i/25/T c q is the BCS coherence length, V 
is the diffusion coefficient. The fluctuating order pa- 
rameter A^-* (r, t) varies close to T c o on a larger scale 
£gl (e) ^> Cbcs- The ratio of the FCP concentration to 
the corresponding effective mass with logarithmic accu- 
racy can be estimated as n c . p ./m c .p. ~ £01° ( e ) and in 
the 2D case assume as canstant (which is the case we will 
discuss in the following)!!!. 

The two principal fluctuation contributions to conduc- 
tivity close to T c o, are positive and originate from a direct 
FCP charge transfer (AL contribution) 



5a 



AL 



(n c .p./m c .p.) e 2 r G L ~ e 2 /he 



(8) 



and from the specific quantum process of one-electron 
charge transfer related to coherent scattering of electrons 
on elastic impurities, which leads to the formation of 
FCPs (anomalous MT contribution) 



j- In (e/ 7 ,). 



However, these two contributions do not capture the 
complete effect of fluctuations on conductivity. The in- 
volvement of quasi-particles in the fluctuation pairing re- 
sults in their absence at the Fermi level, i.e., in the open- 
ing of a pseudo-gap in the one-electron spectrum and 
consequently decrease the one-particle Drude-like con- 
ductivity. Such an indirect effect of the FCP forma- 
tion is usually referred as the DOS contribution. Be- 
ing proportional to the concentration of the FCPs n c p ., 
the DOS contribution formally appears by integration 



of the Fourier-component /|AW (q, w)| \ of the order 
parameter over all long-wave-length fluctuation modes 
(q < CbcsV^); m the static approximation (lj —t 0) given 
by: 



5a 



DOS 



2n r 



e 2 T 



e 1 1 

— In-. 

n e 



(9) 

One sees that the DOS contribution has the opposite sign 
with respect to the AL and MT contributions, but close 
to T c o does not compete with those, since it turns out to 
be less singular as a function of temperature. 

Finally, the one-electron diffusion coefficient is renor- 
malized in the presence of fluctuation pairing (DCR). 
Close to T c q this contribution is not singular in e (see 
table |) and was usually ignored in literature, but as was 
mentioned before, it becomes of primary importance rel- 
atively far from T c q, and at very low temperatures. It 
is due to 5a^ K that the sign of the total contribution 

of fluctuations to conductivity Sai^ changes in a wide 
domain of the phase diagram and in particular close to 
T = 0, in the region of quantum fluctuations (see Fig. |^, 
where the regions with dominating fluctuation contribu- 
tions to magneto-conductivity are shown). 

At zero temperature and fields above H C 2 (0), the sys- 
tematics of the fluctuation contributions to the conduc- 
tivity changes considerably with respect to that close to 
T c q. Due to the collision- less rotation of FCPs (they do 
not "feel" the presence of elastic impurities, all informa- 
tion concerning electron scattering is already included 
in the effective mass of the Cooper pairs) they do not 
contribute directly to the longitudinal (along the applied 
electric field) electric transport (analogously to the sup- 
pression of the one-electron conductivity in strong mag- 
netic fields (lu c t >• I): 8a x °x ~ (lu c t)~ 2 , see Ref. [pl| ) 

and the AL contribution to Saxx ^ becomes zero. The 
anomalous MT and DOS contributions tend to zero as 
well but because of different reasons. Namely, the former 
vanishes since magnetic fields as large as H C 2 (0) com- 
pletely destroy the phase coherence, whereas the latter 
disappears since magnetic field suppresses the fluctuation 
gap in the one-electron spectrum. Therefore the effect of 
fluctuations on the conductivity at zero temperature is 
reduced to the renormalization of the one-electron dif- 
fusion coefficient. FCPs in the quantum region occupy 
the lowest Landau level, but all dynamic fluctuations in 
the frequency interval from to Abcs have to be taken 
into account. The corresponding fluctuation propagator 
at ze ro te mperature close to H C 2 (0) has the form (see 
Eq. (|A22|) ) 



L (oj) = -u 1 



1 



h + cj/Abcs 



and 
5a 



DCR 

xx 



Abcs 



Abcs 



die 



A E 



e 1 I 

--ln T . (10) 
ft h 
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The parameter h = [H — H C 2 (0)] /H C 2 (0) plays the same 
role as the reduced temperature e in the case of the clas- 
sical transition; Abcs is the BCS value of the gap at zero 
temperature in zero field. 

While the denominator of the integrand in Eq. (^) 
defines the characteristic wavelength £gl (T) of the fluc- 
tuation modes close to T c o, the one in Eq. @ defines 

the characteristic coherence time tqf (h\ of QFs near 

H C 2 (0) (where t -C h). The value of the integral is de- 
termined by its lower cut-off wqf ~ Abcs' 1 : an d the 
corresponding time scale is 



(Abcs/i) 



(11) 



One sees that the functional form of tqf is completely 
analogous to that of tql: Abcs — T c0 and the reduced 
field h plays the role of reduced temperature e. Eq. 
( pT| ) can also be obtained from the uncertainty princi- 
ple. Indeed, the energy, characterizing the proximity 
to the quantum phase transition is AE = ftiu c (H) — 
{H C 2 (0)) ~ Abcs^- and namely this value should be 
used in the Heisenberg relation instead of ks(T — T c q), 
as was done in the vicinity of T c q. The spatial coher- 
ence scale £qf (jij can be estimated from the value of 

tqf analogously to the consideration near T c q. Namely, 
two electrons with coherent phase starting from the same 
point get separated by the distance 



C QF (hj ~ (Dt qf ) 



1/2 



i 



BCS 



after time tqf- 

To clarify the physical meaning of tqf and £qf, note 
that near the quantum phase transition at zero temper- 
ature, where H — > H C 2 (0), the fluctuations of the order 
parameter A^ fl ' (r, t) become highly inhomogeneous, con- 
trary to the situation near T c q. Indeed, below H C 2 (0), the 
spatial distribution of the order parameter at finite mag- 
netic field reflects the appearance of Abrikosov vortices 
with average spacing [close to H c2 (0) but in the region 
where the notion of vortices is still adequate] equal to 



a (H) = £bcs /y/H/H& (0) -> £bcs- 

Therefore, one expects that close to and above H c2 (0) 
the fluctuation order parameter A^ (r, t) also has a 
"vortex- like" spatial structure and varies over the scale 
£bcs an d being preserved over time tqf- In the lan- 
guage of FCPs, one describes this situation in the follow- 
ing way: A FCP at zero temperature and in magnetic 
field close to H C 2 (0) rotates with Larmor radius tl ~ 
wfK (H c2 (0)) ~ wf/Abcs ~ £bcs, which represents 
its effective size. During time tqf two initially selected 
electrons participate in multiple fluctuating Cooper pair- 
ings maintaining their coherence. The coherence length 

£qf [hj ^ £bcs is thus a characteristic size of a cluster 
of such coherently rotating FCP, and tqf estimates the 



lifetime of such a flickering cluster. One can view the 
whole system as an ensemble of flickering domains of co- 
herently rotating FCP, precursors of vortices (see Fig. ||) . 

In view of the qualitative picture of SFs in the regime of 
the QPT, let us continue with the scenario of Abrikosov 
lattice defragmentation: Approaching H c2 (0) from be- 
low, puddles of fluctuating vortices are formed, which 
are nothing else as FCPs rotating in a magnetic field. 

Their characteristic size is £qf , and they flicker 

in the characteristic time tqf M/i|V In this situation, 

the supercurrent can still flow through the sample until 
these puddles do not break the last percolating super- 
conductive channel. The corresponding field determines 
the value of the by QFs renormalized second critical field: 
H* 2 (0) = H c2 (0) [1 - 2Giln (1/Gi)] (see Ref. §). Above 
this field no supercurrent can flow through the sample 
anymore, i.e., the system is in the normal state. Nev- 
ertheless, as demonstrated by the above estimates, its 
properties are strongly affected by the QF. Fragments of 
the Abrikosov lattice can be still observed in this region 
by the following Gedanken experiment: The clusters of 
rotating FCP ("ex- vortices") of size £qf with some kind 
of the superconducting order should be found in the back- 
ground of the normal state, if one takes a picture with ex- 
posure time shorter than tqf. For exposure times longer 
than tqf, the picture is smeared out and no traces of 
the Abrikosov vortex state can be found. However, the 
detailed nature of the order which exists there is still un- 
clear. It would be attractive to identify these clusters 
with fragments of the Abrikosov lattice, but most proba- 
ble this is some kind of quantum FCP liquid. Indeed, the 
presence of structural disorder can result in the forma- 
tion of a hexatic phase close to H* 2 (0), where the trans- 
lational invariance no longer exists, while at the same 
time conserving the orientational order or the vortices. 



VI. DISCUSSION 

In terms of the introduced QF characteristics tqf and 
£qf, one can understand the meaning of already found 
microscopic QF contributions to different physical values 
in the vicinity of H C 2 (0) and derive others which are 
related. 



A. In-plane conductivity 

For example, the physical meaning of Eq. (Q) can be 
understood as follows: one could estimate the FCP con- 
ductivity by merely replacing tgl — > tqf in the classi- 
cal AL expression (||), which would give 6a AL ~ e 2 TQF- 
Nevertheless, as we already noticed, a FCP at zero tem- 
perature cannot drift along the electric field but only ro- 
tates around a fixed center. As temperature deviates 
from zero, FCPs can change their state due to the in- 
teraction with the thermal bath, i.e. their hopping to 
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an adjacent rotation trajectory along the applied electric 
field becomes possible. This means that FCP can partic- 
ipate in longitudinal charge transfer now. This process 
can be mappecL-Qnto the paraconductivity of a granular 
superconductorca at temperatures above T c q, where the 
FCP tunneling between grains occurs in two steps: first 
one electron jumps, then the second follows. The proba- 
bility of each hopping event is proportional to the inter- 
grain tunneling rate T. To conserve the superconducting 
coherence between both events, the latter should occur 
during the FCP lifetime r GL . The probability of FCPs 
tunneling between two grains is determined by the con- 
ditional probability of two one-electron hopping events 
and is proportional to Wr = T 2 tql- Coming back to 
the situation of FCPs above H c2 (0), one can identify 
the tunneling rate with temperature T while tql cor- 
responds to tqf- Therefore, in order to obtain a final 
expression, 5a Ah should be multiplied by the probability 
factor Wqf = t 2 TQp of the FCP hopping to the neigh- 
boring trajectory: 



Sa^ ~ 6a AL W QF ~ e 2 t 2 /h 2 , 
which corresponds to the asymptotic Eq. (|]). 

B. Magnetic susceptibility 

In order to estimate the contribution of QFs to the fluc- 
tuation induced magnetic susceptibility of the SC in the 
vicinity of H c2 (0), one can apply the Langevin formula 
to a coherent cluster of FCPs and identify its average size 
by the rotator radius. One finds 



AL 



£bcs/ ch- 



in complete agreement with the result of Rcf. |TE 



C. Nernst coefficient 

One further reproduces the contribution of QFs to 
the Nernst coefficient. Close to H c2 (0) the chemi- 
cal potential of FCPs can be identified as /ifcp = 
huj c (H c2 (0)) - hw c (H) [as in Rcf. [jTJ|, close to T c0 , 
Mfcp = &b {Tco — T)] . The corresponding derivative is 
dfipcp/dT ~ dH c2 (T) jdT T/A BC s- Using the rela- 
tion between the latter and the Nernst coefficient, it is 
possible to reproduce one of the results of Ref. (ll| : 



,al 



[tqfA 



ic.p.] dfiFCp/dT ~ £ BCS 



t/h. 



D. Transversal magneto-resistance above H C 2 (0) 

The proposed qualitative approach can also explain 
the non-monotonic behavior of the transversal magneto- 



resistance observed in the layered organic superconduc- 
tor k — (BEDT p TTF) 2 Cu(NCS) 2 above ff c2 (0) at 
low temperatures^. Indeed, the motion of FCPs along 
the z-axis in such a system has hopping character and the 
quasi-particle spectrum can be assumed to have the form 
of a corrugated cylinder. Close to T c o the fluctuation 
magneto-conductivity tensor in this model was already 
studied in details in Ref. || . There it was demonstrated 
that the transverse paraconductivity in that case is sup- 
pressed by the square of the small anisotropy parameter 
(£,z/£,x) 2 , while the dependence on the reduced tempera- 
ture e is even more singular than in plane. In terms of the 
Ginzburg-Landau FCP life-time (^) , it can be written as 



2 C 4 



^zz (6) 



(12) 



where s is the interlaycr distance. In principle this result 
could be obtained, even from the Drude formula applied 
to the FCP charge transfer [see above, how Eq. (§j) for 
5<?xx ( e ) was obtained] combined with the above specula- 
tions regarding the hopping of FCPs along z-axisc3. This 
general approach, which does not involve the GL scheme, 
allows us to map Eq. ( |l2| ) on the case of the QPT by just 

replacing t gl (e) -> r QF (h) 



5a 



AL 



CO 



4 e2 £z rp 2 

£2 S 3 J cO r QF 



V / £2 S 3 V 7T / h 2 



The negative contribution appearing from the diffu- 
sion coefficient renormalization competes with the posi- 
tive <5er^ z L Chj . The only difference between the in-plane 

[see Eqs. (0) & (10)] and z-axis components of this one- 
particle contribution consists in the anysotropy factor 
( v z) l v x = QJ^i- As a result one gets: 



DCR) 

zz 



2 = 2 «M„i 

h 



and the total fluctuation correction to the z-axis magne- 
toconductivity at zero temperature above H c2 (0) can be 
written as 



2 £2 



2e 2 C 
3tt 2 £ 



1.94 



1 



In- 



1 



(13) 



We used Eq. (13) fat-, the analysis of unpublished 
data by M. KartzovnikE3 on the magneto-resistance 
of the layered organic superconductor k — (BEDT — 
TTF) 2 Cu(NCS) 2 at low temperatures and magnetic 
fields above H c2 (0). The measurement was taken at 
T = 1.7K with a T c0 « 9.5K and B c2 (0) sa 1.57T and 

this curve was fitted by 0.23 (o.l8/h 2 + hihj , see Fig.0. 

For the material parameters of this compound, the au- 
thor reports r = 1.7ps, £ z = 0.3 — 0.4 nm, and s = 1 
nm. The fitting shown in Fig. [l^ corresponds to the ra- 
tio £z/s = 0.32 and looks rather convincing. 
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—■ measurement, T=1.7K 
— fit, t=0 with B (0)=1.57T 
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FIG. 10. (Color online) Comparison to resistivity measure- 
ments of the layered organic superconductor k — (BEDT — 
TTF) 2 Cu(NCS) 2 |||. The material has a transition tem- 
perature of T c0 w 9.5K, B c2 (0) ~ 1.57T, and r = 1.7ps. This 
experimental curve is taken at T = 1.7K and fitted by ex- 
pression in Eq. dli|), which is in perfect agreement with the 
experiment. Specifics are given in the text. 



The discrepancy appearing between the theoretical 
and experimental curves in the high field region, M. 
Kartsovnik attributes to the large normal-state magneto 
resistance, reflecting the specifics of the cyclotron orbits 
on the multi-connected Fermi surface of the compound 
(due to the low crystal symmetry it is quite difficult to 
fit). 
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expression is 



O fc {n,m}=0 

(Al) 

x B%l(ci k ,n k+l/ )L n (n k+v ). 

The block of three Green functions B nm with velocity 
operator (originating from the current vertex) and two 
Cooperons is given by 

B nm (n k+V , ft*) = T^T Tr {G (e,:) (e i+u ) x 

A n (ei+i/, Cl k -i)G A m (Ofc_i, £j) j . 

(A2) 

The trace operator Tr denotes the integration over all 
electron quantum numbers. The corresponding block was 
calculated in |Q exactly for fields with lu c t <C 1, i.e. for 
the case of our interest. Under this condition the Landau 
quantization affects the motion of Cooper pairs, while 
the Green functions in the block Eq. (A2) can be used in 
r— approximation. As the result, using the properties of 
the velocity operator in Landau representation, one finds 



fl<>> (fin-*, 0*) = -2u V y/eH(n + l)S m , n +i 

(n k+ „,n k ), (A3) 



+V eHndn 



with 



6 (s i+u ttk-i) 



\2ei+w v -n k \+iu c (n + 1/2) 



n k \+uj c {m. + l/2)' 



Appendix A: Aslamazov-Larkin contribution 

1. General expression 

Let us start with the discussion of the AL contribu- 
tion (diagram 1 in Fig. [j]). The corresponding analytic 



Substituting Eq. (A3) in Eq. (Al) and further summa- 
tion over Landau levels in Eq. ( Al ), results in the cancel- 
lation of the terms containing the products <S m ,n+i<5n,m+i 
and S m . n iS n.m-i- The analysis of the theta-functions 
in Eq. (A4) results in the possibility of separation of 
different domains of analyticity in the plane of bosonic 
frequencies tt k : 
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FIG. 11. The integration contour in the plane of complex 
frequencies. 



S m „(ri fe , fl k + U„) =27tT 



fe-l 



i— — OO i — — DC 

1 



i—k 

l 



i=0 



|2e< +w„ - fi fc | +w c (n+ 1/2) |2ei - f2 fe | +w c (m+l/2) 
Summation over fermionic frequency in this expression can already be performed in terms of ip— functions: 

1 



1 ^ + |» fc |+^ c (n + l/2) \ _ /l |n fc |+u; c (m+l/2) 
1 2 AttT J ^ V 2 4ttT 



2w c (n — to) 

'1 |fi fc+I/ | +w c (n + 1/2) 



1> 



4ttT 



, .'1 ^ + |n fe+y |+a; c (m + l/2) 
1 2 4ttT 



(A5) 



(A6) 



Being interested in the d.c. fluctuation conductivity, i.e. 
taking into account the li mit uj v — > — iui — > after analyt- 
ical continuation, in Eq. (|A6| ) we neglected the frequency 
w„ in comparison with uj c (n — to) in denominator since 
the diagonal term (to = n) disappears in the process of 
summatio n ov er Landau levels in Eq. (Al) as follows 
from Eq. (A3). One notices the useful fact that the per- 
mutation fifc fife + uj u simultaneously with to n in 
Eq. (|A6|) does not change the function S mT[ (Qfc, f2fc+w, y ) : 



(A7) 



Let us retur n t o the general expression for paracon- 
ductivity Eq. (Al). One can transform the sum over the 



bosonic frequencies flk to the contour integral / in the 
plane of complex frequency i7fe — > —iz: 



QxxM = -We^V'eH V C mn I^ n («„) , (A8) 



AL 



C K) = — j cotn (^;) dzZ nm (-iz + u v ,-iz)x 

(A9) 

3ron(-K, — + 0J v )L m {—iz)L n (-iz + W„), 

where the contour integral encloses all frequencies fi^ [in 
the plane of frequency z these are poles of coth(z/2T), 
see Fig. |ll). The coefficients 



(^m,n+l^n,m — 1 ~\~ ^n,m-\-l^m y n— 1 

(A10) 

control the summation over Landau levels. 



Let us stress that both functions S in Eq. (A9) have 
breaks of their analyticity along the lines Imz = and 
Imzz = the same as the product of the propagators. 
As a result, one gets three domains where the integrand 
function is analytical: above the line Imz = 0, between 
the lines Imz = and Imzz = — w v and below Im z = 
— uj v . For the analytical continuation of function ( |A6| ) 
to the whole complex plane from Matsubara frequencies, 
three different functions: S^,S^, and S^, should be 
introduced, which are analytical in their corresponding 
domains. They differ by the combinations of the signs of 
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the explicit absolute values appearing in Eq. (A6). Due 
to observation ( A7 ) one can write the useful identities 



-rr 



a^Hz + w„, -iz) = E%£(-iz, -iz + u v ) 



RAi 



and get for the contour integral in Eq. (A9) 



Aml^ K) = J coth (^) dz { [E™(-iz + u v , -iz)} 2 L*(-iz) - [Z%£(-iz + u v , -iz)} 2 L^(-iz)} L*(-iz + u v ) 

—iw,, / Z \ ( 2 2 1 

coth (-)dz{[Z™(-iz + uj u ,-iz)] LK(-iz + u v )-\p™{-iz + Lo v ,-iz)\ L*(-iz + uj u )j L*(-iz). 



The last integration can be r educed to that along the real axis by means of shifting the variable — iz + uj v — > —iz' . 
The resulting expression (A9) for the electromagnetic response operator - still defined on Matsubara frequencies uj v 
- takes the form: 



t,2 eH c ^ 

7T ^ " 



Q^(uj v )=^e 2 u 2 V 2 



coth 



$ m „ (z,u v ) dz, 



(All) 



where 



4>, 



(z, u v ) = { [S*£(-w + w„, -«)] 2 (-»«) - [S££(-« + u u , -iz)] 2 Lf n (-iz)} L R (-iz + uj v ) 

+ { [ZmiHz - w„, -iz)] 2 L*(-iz) - [s££(-iz, -i* - w„)] 2 L*{-iz - «„). (A12) 



The rules for performing the analytical continuations of the function H mn (Ofe, Clj. + in Eq. ( A12 ) are simple: the 
sign of the explicitly written absolute values of the corresponding frequency in Eq. (A6) is chosen as "+" in the case 
of retarded continuation (superscript R) and it is chosen as "— " in the case of the advanced one (superscript A). For 
instance 



2uj c (n — m) 



1 uj v - tt k +uj c (n + 1/2) 

2 4ttT 



-n k +uj c {m + l/2) 
4ttT 



'1 w„+n fc + w c (n+l/2)\ . /I , 2w v + n fc + w c (m+ 1/2) 



4ttT 



4ttT 



and analogously for and 

Now one can perform the last analytical continuation 
uj v — > —iuj m Eq. ( |A12|) and obtain (z.w) as an 
analytic function of the real external frequency uj. Since 
we are interested in the d.c. limit of the FC, i.e. u — > 0, 
the function <E>^?2 ( z i w ) can be presented in the form of 
its Taylor expansion: 



to 2 (n — rn) 1 



-iz). 



The first term is not of interest here: all frequency inde- 
pendent contributions which form (0, T, H) are can- 
celled out: this is a necessary requirement of the absence 
of the diamagnetic response in the normal phase of super- 
conductors. Actually, in order to find the FC we need to 
know only ImQ^ (uj, T, H) , i.e. we are interested only 



I 

the imaginary part of F „ m (-iz). It can be obtained by 
expansion of all functions H"^ (a, (3 = R,A) and propa- 
gators in (z,uj) over uj. Introducing the function 



^nm(iz) = ip 
-1> 



1 iz + uj c (n + 1/2) 

2 + Arff 

1 iz + uj c (m+ 1/2) 

2 + 4^T 
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one can find the an alytic ally continued expressions for 
the products of Eq. ( |A12| ): 



L^nmJ 



[E*£(-iz±iu, 
_ ReV nm (-iz) 
w 2 (n — to) 2 
\^AA\ 2 

L 1 — 'nmj 



w c ( n ~ m ) 

■iz)] 2 

Re* nm (-iz) ± -^^' nm (±i z ; 



1 



w 2 (n — m,y 



4ttT 
(iz) + O (w 2 )] 



which leads to 
ImF „ ro (-iz) 



<9z 



{2Re* 2 lTO ImL«Imi« 



+Im* 2 m [lmZ*Re£* + ImL*ReL«] } . 

One can see that this function is symmetric with 
respect to subscripts permutation: ImF nm (— iz) = 
ImF mn (- iz). Let us stress that we have could present 

the linear in oj part of the function (z, oS) in the form 
of full derivative with respect to z. The same situation, 
was found in the original paper of Aslamazov and Larkinel 
for the simple case when the Green functions block could 



be assumed to be constant. As a consequence of this im- 
port ant p roperty of <&mn (z,uj), the integration over z in 



Eq. (All) can be perform ed b y parts. After summation 
over to, the Eqs. (Q) and (AS) read as 



8<r£{T,H) 



vZYXn + 1) 



2-kT ^ 

n=0 



dz 



sinh 2 (z/2T) 
{2Re* 2 lj „ +1 (-^)ImL«(-w)Imi« +1 (-w) 



Im^,„+i(-^) [ImL%(-iz)ReL% +1 (-iz) 
ImL% +1 (-iz)ReL%(-iz)]} . 



(A13) 



Let us attract the attention to the fact that due to the 
integratio n by parts coth z/2T disapp eared from the inte- 
gral Eq. ( All ) being replaced in Eq. ( |A13| ) by its deriva- 
tive sinh - (z/2T) . This fact makes our answer different 
from the one of Ref. |l5| and physically means, as we will 
see below, that at low temperatures the paraconducting 
contribution tend to zero: fluctuation Cooper pairs above 
H C 2 (0) exist but do not moye-and do not participate di- 
rectly in the charge transfcrtHiil. 

It is convenient to introduce the dimcnsionless variable: 
x = z/(2nT), parameters t = T/T c0 and h = 2e£ 2 H, 
where £ 2 = irT)/ (8T) , and the function 



1 + ix 2 fh\ , 



£ m (t, h, ix) =\nt + ip 
In this representation, Eq. ( |A13| ) takes the form (£k (t, h,ix) = £k) 

<fo£ (*,*) = ^£> + l)/ 

™ — n J — 



'5 



'_no sinh TTX 



dx j Im 2 £ m ^ Im 2 g m+ i Im 2 £ m +i - Im 2 £ m t> \c c 1 

\£m\ \£m+l\ 



\£. 



(A14) 



(A15) 



m+l| 



This is the general expression for fluctuation paracon- 
ductivity valid in all domains of temperatures and fields 
under consideration. 

We will see that all diagrams presented in Fig]l] are rel- 
evant in different regions of the phase diagram, depicted 
m Fig. |. Nine regions of different asymptotic behav- 
ior can be distinguished and below we will analyze all 
contributions in each domain. 



2. Asymptotic behavior 

a. Vicinity of T c0 , fields h < 1(H < H c2 (0)) 

I n thi s case hit = e <C 1 and the -^-function in 
Eq.(Al4) can be expanded. In first approximation: 



£W {t ,h,ix) = e+ l ^+(^)(m+l). (A16) 



The integral in Eq. (A15) can be easily carried our: only 
the first fraction in the parenthesis should be taken into 
account. Further summation over Landau levels can be 
performed exactly in terms of the ^-function: 



^M<<1) = ^UJ 



1 l + sk 



(A17) 

which coincides with the known expression for the 
Cooper pairs contribution to the magneto-conductivity 
in the Ginzburg-Landau region!!. 

The general Eq. (A15) allows to obtain the next order 
correction in e with respect to the AL result. In order 
to do this, one should take into account both terms and 
expand up to the second order: 

£%> (t, h, fa) = £&) (t, h, ix) - (to + \ 

m+l) (A18) 
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After some simple but cumbersome calculation in the 
limit of small fields, one finds 



^i L (e«l) 



e 

16e 



7C(3)e 2 1 1 



(A19) 



which demonstrates that the next order correction to the 
well known AL result in the vicinity of T c o is smaller than 
the DOSu and the regular MTEJ contributions (see below) 
by only numerical prc-factor. 



b. High temperatures T S> Tco , weak fields h <C t 

Let us move to the discussion of the high-te mper ature 
asymptotic. We will assume hit ^> 1 in Eq. (A14) and 
get: 



"1 






2- 







Im£ TO (t, h, ix) = ~ip' 

J (A20) 

The sum in Eq. (A15) converges at n max ~ t/h 1 
and can be replaced by an integral. The integration over 
x involves only the region x ~ 1 and can be performed 
first. As a result one gets 



8<r£x(t > l,h < t) 



(>7T 



ln 3 i 



with Ci = | J °° [V'' (1/2 + x)] 3 dx = 6.97. Let us stress 
that this asymptotic expression coincides with the high 
temperature behavior of the AL contribution obtained in 
clean caseE3 which emphasises the statement that the 2D 
paraconductivity is an universal function of In t through- 
out the complete temperature range. 



c. Fields close to the line H C 2 (T) 

The line separating normal and supcrconductiving 
phases H C 2 (T) (in our dimcnsionlcss units the line of 
critical fields h C 2 (t)) is determined by the requirement 
that the propagator (||) has a pole when Slfc = and 
m = 0: 



In t + if) 



2 h c2 (t) 
ir 2 t 



-1> 



= 0. 



At low temperatures T <C T c q, close to the point 
T = and H = H c2 (0), the critical field is h c2 (t) = 
2£, 2 H C 2(0)/e^ 1. Then one can substitute the -0-function 
by its asymptotic expression ip( x ) = In a; — l/(2x) and 
take into account that ip(l/2) = — \nijE {lE = 1.781.. 
is the Euler's constant) which results in 



h c2 (t -> 0) 



87i 



(A21) 



In order to find the paraconducting contribution to 
FC above the curve H C 2 (T) in Fig. let us rewrite Eq. 



( A14) in terms of the reduced field 



Below we will see that the Cooper pair contribution to 
FC, which is singular in ft, -1 , originates in Eq. ( |A15| ) only 
from the term with m = 0, i.e. we can restrict ourselves 
to the Lowest Landau Level (LLL) approximation. Hence 

we will need the explicit expression for £ m [t, h, ix^j only 

for m — 0, 1 and h <t« /i C 2 (t) ■ In order to get this, 
one can use in Eq ( A14 ) a paramctrization in terms of 
h and expand it h <C 1 and h — h C 2 (t) = h • h C 2 (t) <C t. 
This gives 

~ m 2 xt 

h+ -r j-r . A22 

4h c2 (t) 



£o 1 1, h, ix 



The substitution of Eq. flA22[ ) to Eq. ( |A15| ) results in 

1 4h c2 (t)h\ 



8<r£(t,h) 



;Jgl 



■kH 



(A23) 



with 



Jgl (r) 



dx 



sinh x x z 



3-2 = 2W/ (r) - - - 2 
(A24) 

first calculated in Ref. ||15[. This formula is valid along 
all the line h C 2 (t) until t ~ h c2 (t). Taking into account 
the asymptotic expressions 



N 1 1 1 
ip (r^oo) = - + _ + —; 

7* Zr vr 
one finds that in this domain 



V>'(r->0) = 1/r 2 . (A25) 



, 3tt 2 /i 2 



, * < h, 



, h c2 (t) h<^t. 



(A26) 



k 4/i c2 (t)/» 

The first line of Eq. ( |A26 ) corresponds to the quan- 
tum fluctuations which arc realized in the limit of low- 
est temperatures t <C h close to H C 2 (0) . One sees that 
the paraconductivity decays here as T 2 . In the temper- 
ature range h C 2 (t) h <C t <C h C 2 (t) the parac onductivity 
is determined by the first line of Eq. ( A26 ). Close to 
H C 2 (0) but for relatively high temperatures t ~ h the 
corresponding expression can be rewritten using the ex- 
plicit expression for h C 2 (0) Eq. (A21): 



8<r£(t,h) 



(A27) 



which perfectly matches to the first line of the Eq. flA26Q . 
Here the transition from quantum to classical fluctua- 
tions takes place. At higher temperatures along the line 
h C 2 (t) one should take into account the temperature de- 
pendence of h C 2 (t): 

-ts=w (A28) 

This expression is valid along the line h C 2 (t) until t <C 
h c2 (t) where Eq. (|A28[) matches to Eq. flAl7|). 
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d. High fields (H 3> H C 2 (0) ), temperatures t <C h 

In this domain we are far from the transition line 
H C 2 (T), [h ^> h C 2 (t))} and the LLL approximation is not 
applicable. Nevertheless, one can substitute the summa- 
tion over Landau le vels b y an integration. Replacing the 
ip— function in Eq. (A14) by a logarithm, one finds that 



„ . , , Ah f 1 
E m (t, h, ix) = In — I m + - 



'5 



in xt 



Ah (2m + 1)' 
(A29) 

One can see that this expression reproduces Eq. ( |A22| ) 
when h h C 2 (t) a nd m = 0. Let us substitute Eq. 



( |A29| ) to Eq. ( A15 ). As we will see below, the sum 
converges at m ~ 1. It is why th e m ain contribution 
comes from the second term of Eq. ( A15 ) , where the sum 
Im£ m Re£ m +i + Im£ m+ iRe£ rn ~ In fMj . As a result we 
get 



192 



Co 



In 



3 2h 



with C-2 = 0.545. Let us recall that this expression is 
valid for arbitrary temperatures small in comparison to 
reduced field t <C h. 



Appendix B: Maki-Thompson contribution 

1. General expression 

Below we will calculate the fluctuation rcnormalization 
of the one-electron contributions to conductivity. It is 
technically convenient to start with the usual expressions 
for the Maki-Thompson and other diagrams from Fig. [l], 
written in momentum representation. Only at the very 
end one should quantize the motion of Cooper pairs in a 
magnetic field in accordance with the rule 

M 



Diagram 2 from Fig. |l| can be written as 



where 



(Bl) 



E^ T (q, n fc , w„) = A(q, e n+v , ft fc _„_„) 

x A(q,£„,rjfc_„)/^, T (q,£„,^ fc ,a; 1 ,) 



(B2) 



and 



!f7 = J 7^^«o;(p)v a! (q-p)G ! (p,£ n+v ) 

x G(p, e n )G(q — p, Ofe — n — h> )G(q- p, flk-n) 



The main q-dependence in (Bl) arises from the propaga- 
tor and vertices A. That is why we can assume q = in 
the Green functions and calculate the electron momen- 
tum integral by changing, as usual, to a £(p) integration: 



rMT t, -1,, 



1 



1 



1 



(B3) 



iO, 



'At. 



1 c„ 



J £ n+ , 
> i £ n 



FIG. 12. ^-integration in Eq. (B3): all six possible positions 
of the poles in the complex plane of £ are shown. 



This integral, Eq. (B3), can be calculated using the 
Cauchy theorem. Closing the contour in upper or lower 
half-plane by the large semicircle and noticing, that, 
due to fas t decrease of the integrand the function in 
Eq. (B3), the integral over the semicircle becomes 
zero, one can express J xx in terms of the sum of 
the corresponding residues. There are six different 
combinations of the pole positions with respect to the 
real axis in the c omp lex plane of £, leading to non-zero 
results (see Fig. [l|): two realization corresponding to 
9 (— £„£„ + „) 9 (Clk-n^k-n-v) 7^ 0, one realization corre- 
sponding to 9 (-£„£„+„) 9 {Sl k - n VL k - n - v ) ^ 0, two re- 
alization corresponding to 9 (£ n £n+v) 9 (Qk-nQk-n-v) 7^ 
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0, and the realization corresponding to 
9 (£„£„+„) 9 (-£l k - n n k - n -v) ^ 0. Calculating the 
residues for each situation and assuming that 
e n = (2t) 1 sgn£„ (let us recall that we consider 
the disordered limit T <C one hnds: 



I™ 1 = 2irVis T 2 {9 (s n e n+v ) 9 (fl k _ n Cl k _ n _„) 

+ [£n £ n+v) 9 {— fifc— n^fe-n-i/) 

— 29 (—e n e n +u) 9 {—^k-n^k-n-v) 

— 29 (£„£„+„) 9 (Qk-n^k-n-v)} ■ 

Now one should substitute this expression to Eq. ( p32| ) 
and perform the summation over the fermionic frequen- 
cies. This is a cumbersome exercise, which, nevertheless, 
can be followed through analytically. Here we mention 
some useful transformations which are important to per- 
form the summations: One can see, that the simultane- 
ous permutations n — > —n and k — > — k allows to simplify 
the sums: 



E£ T = (e£ t (°») + S^ T(rcg2) ) + £™ T(regl) = -2w VT 



^ n=— i 

CO 



29 (— rife_ n Ofc_„_y) — 9 (rt k - n n k ^ n - v ) 

\e n +u - Slk-n-v\ +Vq 2 ) (|e„ - f2fc_„| +Vq 2 ) 

29 (Ofc_„Ofc_ ra _j,) — 9 ( — 17fc_„r2fc_ Jl _ l/ ) I 

(\s n+v - n fc _ n _„| + v q 2 ) (| e „ - n k _ n \ + 2V) f 



The rules writing the absolute values explicitly in the 
sum using the first 9- function is evident. In the sec- 
ond sum, containing 9 (p, k - n VL k - n ^ v ) , one should make 
a shift n' = n + v . After rewriting the absolute values 
for the Cooperons, the sums can be expressed in terms 
of ^-functions. Using the identity 



V(l/2- 



^(1/2 



ni tanh ttz 



allows to write the final expression for the first sum as 

yMT(an) , yMT(rcg2) _ P^O^ {u v -i - \ fl k \) 

, / 1 2u„- \Q k \+Vg 2 \ _ (\ \n k \+v q 2 ^ 
^\2 4ttT J ^\2 4ttT 

(B4) 

Looking at the denominator of this expression one can 
recognize that E^ c T ^ an ' ) is responsible for the anomalous 
Maki-Thompson term. 

Next we consider the remaining second sum in S^ T . 
One can see that in the first term both arguments of 
the absolute values are positive. In the second term 
we can replace k — > —k, with an additional change of 
the order of the summation over bosonic frequencies. 
The sum with 9 {V Lk-rSlk-n-v) can be calculated in the 
spirit of Eq. (B4). Regarding the last sum, containing 
9 {—VL k ^ n Qk-n~v) , one can find that it is exactly equals 



to zero for any Q k . Finally 



yMT(regl) 







1 , 2w„ + [Ofel +Vq 2 



2 4ttT 

1 \n k \+v q 2 ' 

2 4ttT 



(B5) 



Using the explicit Eqs. (B4)-(B5) we can perform the 



final summation over bosonic frequencies in Eq. (Bl) 
and the analytical continuation of Q^J (w„) to the 
axis of real frequencies. The analytical continuation of 

Q-ra esl ^ is trivial since Eq. (B5) is the analytical func- 
tion of uj v . As a result we get 



)MT(rcgl)R 

'XX 



d 2 q 



J 4tt 2 T J (2tt) 2 

oo , 



h=—oo 



\n k \+Vq 2 
4irT 



Next we go over from the integration over the momentum 
of the Cooper pair center of mass q to the summation 
over Landau levels. Recalling that the density of states 
at the Landau level is H/$q, one finds 



5(T MT(regl) 



where Em ^ (t, h, z) = 9f£ m (t, h, z), i.e 



e 2 VH 



4£^(t,h,\k\) 
^ T $ o f ^ ^(t,h,\k\)" 



E E 



(B6) 



£^(t,h,\k\)=^ 



l + \k\ 2h 



2 + ^t {2m+1) 



In the part of the electromagnetic operator related to 
Eq. ( B4 ) , the external frequency w„ appears in the upper 
limit of the bosonic sum: 

n MT(an) , n MT(reg2) 
'till ' ^xx 

-2e 2 TVv, I f_l ^ 9 J2 L (^k) 



(2tt) 2 w„ + Vq 2 



fc|=0 



i 2^-\n k \ + v q 2 

1 2 AttT 



1 2 4ttT 



and the procedure of analytical continuation is more so- 
phisticated. First of all one can easily see that the contri- 
butions of positive and negative k are equal. The method 
to continue such a sum the real frequencies was developed 
in Ref. j?j and consists in an Eliashberg transformation 
of the sum over tt k to an integral over the contour C (see 
Fig. and the detailed description of this procedure in 
Ref. if). One finds 
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Im z 




Rez 



FIG. 13. Integration contour used in the analytic continuation 
of the MT contribution. 



lMT(an) 
' xx 



Q MT(rcg2)R = _^ TVvq 



+ — i> dz coth (ttz) L(q, — iz) 
2i ' 
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2uj w + Vcf 
4?rT 



c 2 

1 

2 







-il, 
1 

2 ' 



2 Z 

~2 



4ttT/ 



1 



2 4ttT y 
2uj u + £>g 2 
4ttT 



X>g 2 



iz 
~2 



The residue at the point z — iv is equal to zero. Shifting 
the variables in the integral over the upper line Imz = v 
as Z\ = z' + iv, one can present as an analytical 

function of cj„ and analytically continue it in the stan- 
dard way iuj v — > u> — > 0. Expanding it in small w and 
integrating by parts, one gets 



c MT(an) , c MT(reg2) 

T.T. > UU T.T. 



— V 

2 



d 2 q 



(2tt) S 



Vq 2 



dz 



— I z 
2 



4ttT 



2 



4ttT 



sinh ttz 



ln£ 



4ttT 



This expression is then transformed to the summations over Landau leve ls and written in the dimensionless variables, 
in the same way as it was done before. Adding the regular part Eq. ( |B6| ), we finally write the most general expression 
for the MT contribution valid in all domains of temperatures and magnetic field under consideration: 



X MT 
OCT 



r MT(regl) , / r_MT(an) , x MT(rcg2) 

UU XX I \ UU XX ' UU XX 



M 



m=0 kfe=— 00 



Im 2 £ m (i, h, iz) 



£ m (t, h, \k\) 1<p + ( m + 1/2) 7_ OC) sinh 2 ttz Re 2 f m (t, fe, iz) + Im 2 £ m (t, h, iz) 



(B7) 



where M = (tT c or) 1 is the cut-off parameter. 



2. Asymptotic behavior 



Contribution Sa^^" 



Let us start with the eva luation of the contribution 
^ T(r ° sl) given by Eq. ©. 

Vicinity o/T c0 , fieldsh<t: 1(H < if c2 (0)). - Here 
one can just use Eq. (A14) for integer argument k, con- 



sidering the smallness of h/t, expand the ^-function, and 
perform the summation exactly: 



<^ T(resl) (e«l,h) 



7C(3)e 2 



1 et 

2 + 2h 
(B8) 



High temperatures T 3> T c o, weak fields h < C t . 
For the high-temperature asymptotic of Eq. (B6), we 
assume In t 1 . The sum over k is determined by 
£" n (t,h, \k\) and converges fast: it can be performed in 
first. The remaining sum over Landau levels slowly di- 
verges at large m max and can be substituted by an inte- 
gral. The double logarithmic divergence of this integral 
at the upper limit should be cut off at the limit corre- 
sponding m max ~ (T c qt) 1 which results in 



X MT(rcgl) = _'_ 

wu XX 9 



In hv 



T c0 t 



- In hit 



(B9) 
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One can see that close to T c q In (In In t) — > In (7), Eqs. 
(B9) and ( |B8| ) therefore match each other. 

Fields close to the line H C 2 (T) . — In this domain, as 
above, one can use the LLL approximation. Nevertheless, 
since t <C h C 2 ft), many terms are involved in the sum 
over k in Eq. (B6). Therefore it is convenient to present 

the function £q It, h, \k\ ) in this case in the form 



Sb(t,ft,l*l) = 



2 h c2 (t) 

7T 2 t 



1 + h (t) 



ir 2 t 



we first perform the summation over k and integrate over 
Landau levels: 



( 5 cr MT(regl 



— f In In ■ 



T c0 t 



In In — 



(B12) 



The only difference between Eqs. (JBJJ) and (B12) consists 
in the lower limit: in the former it is determined by the 
temperature while in the latte r its role is taken by the 
zero Landau level w c 3> T. Eq. ( B12 ) is valid for arbitrary 
t emp eratures smaller u c and it obviously matches Eq. 
(Bll) along the axis of the magnetic field (t = 0). 



and taking into account the asymptotic expression ( A25) 
one finds 



^MTtrcgl) = _'_ 



2 h(t) \2h 



t 



2e*_ 

7T 2 



In 



7r 2 t T \~ ■ n 2 t 
Along the line H C 2 (T) , in the region of classic fluctu- 



ations h < t -C h C 2 (t) Eq. ( P36|) contains two contri- 
butions: one, originating from the term with k = 0, is 
singular in h (t) — h c 2 (t), the second is logarithmic in 
reduced temperature: 



£ a MT(regl) = _] 



t 



2V ln 2^_ 

4 h- h c2 (t) 7T 2 7T 2 t 

In the the regime of quantum fluctuations f < h, we get 
27_Be 2 t 2e 2 , 1 



^MTXregl) 



■ln^. 

^ h /l 



ffig/i /ie/(ts H 3> i? C 2 (0) , temperatures t -c /i. 
This domain is analogous to the previous one. As above, 



J 



6. Contribution 5(T^; T ' an ' + <5o-^ T ' rcs2 ' 

Now we consider the second part of the MT contribu- 
tion Eq. (B7), namely Sa^^ + 5o-^J^ c&2 \ 

Vicinity of T c0 , fields h < 1(H < H c2 (0)). - In 
the vicinity of the critical temp er atur e T c q one should 
use the expansion Eqs. ( A16 )-( [A18| ) of £ m (t,h,ix). 
For the second order correction it is sufficient to take 
on only the imaginary part of £ m (t, h, ix) into ac- 
count. Substituting correspondingly Rc£m^ (e, ft<l, ix) 
and Imfm' (e, /i <C 1, ix) to Eq. ( p7| ) and using the fact 
that 70 -C 1, the integral over x can be easily performed 
[it converges for x > xo ~ e + (4) (2m + 1)]. The re- 
maining summation is accomplished in terms of the tp- 
functions and its result consists of two terms: the first 
one corresponds to the anomalous MT term So^^^, 



(Bll) while the second, 8a 



MT(rog2) 



gion with 8o~. 
the total So-^J 



MT(rogl) 



, exactly coincides in this re- 



xx - So- xx 
which takes the form 



MT(regi) 



[Eq. (B10)]. Therefore, we present 



do - x 



MT(a 



do x 



MT(rog2) 



p 2 1 
0O- x x = — 



14C(3), 



[ l)-n i/2+ I 



(B13) 



This formula is valid in the vicinity of the critical tem- 
perature T c0 , where we have three different regimes [weak 
fields h <C e, GL strong fields e -C /i, and fields close to 
the /i C 2 (e) line which is "mirrored" at T c q\. 



High temperatures T 3> T c o, weafc fields h <C t- 
Here we discuss the high-temperature asymptotic. As i t 
was done before, we assume hit 3> 1 and use Eqs. ( A14 ) 
and ( A20). Integration over x due to the factor cosh - 7rz 
involves only the region z ~ 1 and can be performed first. 
The sum over Landau levels in this case converges at large 
m max ~f//i> 1 and it can be substituted by an integral. 
The contributing part of the integration with logarithmic 



accuracy turns out to be only the fraction containing 7^. 
As result we get: 



192 ln 2 t ' 



(B14) 



Despite the presence of the large logarithm In in this 
result in its numerator is relatively small with respect to 
Eq. (BS) due to the large ln 2 t in denominator of Eq. 

(|il). 

Fields close to the line H C 2 (T). — As it was done in 



the case of the paraconductivity, let us use in Eq. (B7) 
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the asymptotic Eq. ( A.22Q and perform the calculations 
in the LLL approximation. One easily finds that the 
r esult in this case is also expressed in terms of the integral 
(HH): 



^MT(reg2) = 



7T 2 1 + tjj/h 



Jgl 



4/IC2 (t)- 



(B15) 



We consider the case of low temperatures t <C /i C 2 (t) > 
j and we can omit it in denominator. 



hence 74, < 1 < 

In result it turns out that in this region Eq. ( JB 1 5| ) ex- 
actly coinc ides with the corresponding AL contrib ution 
Eq. ( |A23| ). Therefore it is determined by Eq. ( |A26| ) 
along the whole line H C 2 (T) , i.e. for t <C h C 2 (t), which 
was a lready analyzed in detail above. Looking at Eqs. 
( |A26j ) one can see that in the region h — h C 2 (t) <C t, 
strong cancellation takes place in the MT contribution 
<5<7^ T and only the logarithmic contribution remains 



5a 



MT 



2e 2 2h c2 (t) 
7r 2 ir 2 t 



In the regime of quantum fluctuations t < h — h C 2 (t) 
the contribution Sa^^^ + 5cr x A x TI ' lcg2 ^ ~ t 2 and the lin- 
ear (in t) part of 8(j^^° g1 ^ are gradually frozen and the 
MT contribution reaches the finite negative value at zero 
temperature 



6*™ = — 



2e 2 1 
7T m - 
h 



7T- 



27ge 2 t 
h (t) 



O 



Mt). 



High fields (H H c2 (T) ). — In this domain we are 
far from the transition line H C 2 (T) and the LLL approx- 
imation i s not applicable. Replacing the ^—function in 
the Eq. ( A14 ) by the logarithm one can use the asymp- 
totic Eq. (A29) and gets 



^MTtan) + ^MTfrc^) = 7^3^ 



2 e 2 t 2 



1 



768 h 2 In 2 ^ 



which is beyond the accuracy of the large contribution 
<5cr^ T( ' ros1 ^ [see Eq. (|B12f) l which, in result, determines 
the value of Sa^fJ- in strong fields. 

Finally all asymptotic expressions for the MT diagram 



are summarized in Table II] 



Appendix C: DOS renormalization: contribution of 
the diagrams 3-6 



1. General expression 



We start with calculation of diagram 4. As above we 
use the intermediate results of Ref. for the diagrams 
and then quantize the motion of the center of mass of the 
Cooper pair in magnetic field. The general expression for 
diagram 4 is given by 



(CI) 

with the integral I xx of the four electron Green functions 
calculated exactly in Ref. (!]] in the same spirit as it was 
demonstrated above: 



j(4) 



d 2 p 



v x G (p, e„) G (p, ) g (p, n k _ n ) 



(2n) 2 

= -2-kvqDt 2 [O (e n e n+v ) (e„e„_ fc ) 
+0 (-£„£„+„) (-e„e n _ fe )] . 



Substituting this expression to Eq. (ClT) one finds 



*E 



(e n + n k ) 



E 



,7^0 I 2e « + °* + Vq 2 Y ^ [2e„ + n k + Vq 2 } 



(C2) 



The first term in this expression does not depend on 
external frequency and the corresponding part of the 
electro-magnetic response operator does not contribute 
to conductivity. In the remaining part Q^x(ui v ) one can 
perform the summation over fcrmionic frequency and ob- 



tain it in the form of a sum of two terms: 
v Ve 2 f d 2 q 



4tt 



Em^d 



(\_ Q k +Vq 2 \ ^ + n k + Vq 2 

\2 4irT J ^ I 2 4ttT 



)( 4 ' 2 ) 



K) 



f De 2 /" d 2 q 



4tt 7 (2tt) 2 

' fi n k + vg 2 \ _ n 2uj„-n k + Vq 2 - 
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MT(an) . MT(rcg2)R 

7 II "T ""M 



'(t) 



e j f 

- x- t r rrr 



-/ 2 + &)] 



T" h-h ( .-,(t) 



192 ln 2 j 



?F [*(*)-*(* + *) 



_ X - 2e ,, In 



2/) (f ) 



■ hi 



14C(3)e 



■(*) 



14C(3)e 



■ In 



— rr - 



192 l^r; 



TABLE III. Asymptotic behavior of the MT contributions in different domains, see also Fig. ^ and table [n| 



The analytical continuation of the first one is trivial and 
it gives the first contribution to the conductivity, which 
in Landau representation takes the form 



2e2 

7T 4 



M 



The analytical continuation of Q xx 2 \dJ u ) is completely 
analogous to that one performed above in the case of the 
anomalous MT part. As a result the total contribution 
of diagrams 3 and 4 can be presented as a sum of two 
very different terms 



4e 2 



M 

E 



A ^ (t, ft, A) 7r f°° dx Im£o (t, ft, ix) Ith£' (t, ft, ix) 

£ m (t, ft, k) + 2 J_ 00 sinh 2 (ttx) Re 2 £ (*, ft, ix) + Im 2 £ (t, ft, ix) 



(C3) 



Next, we discuss diagram 5. Its contribution can be 
written in the same way as above: 



)(5) 



K) 



e 2 r 2 i^ /■ d 2 g 

27TI/OT 



/ (^E L ^^ fc )A 2 (?,£„, n fc _ n ) 



where the integral 

d 2 P 



7 (5)( 



(27T) 2 
.2 



u) 



2iriv T sgne n+! ,0 ( ). (C4) 



As result: 

Q^ x (u v ) = -^v Q Ve 2 T 2 



d 2 q 



E i(9,o fc ) 



k— — c 



E 



6 (gfc - En) 
|2e„-fi fe |+2?q 2 ] 



2 ■ 



(C5) 



Further evaluation of this expression is very similar to 
that one of Qxx K). In particular, after summation over 
fcrmionic frequencies, Q xx (lo v ) is presented in the form 



I 

of two sums over bosonic frequencies, one in the limits 
k € [0, oo), the other fee [1, f — 1] and following step-by- 
step the same procedure of the analytical continuation as 
before, one finds that 5a x 5 x ^ = —Scr-kc , and Sai 5 ^ = 
Sa x 4 x 2 \ Therefore we get 



v 7 m=0 



fc=0 ^ m ^' ^' ^ 



(C6) 



7T 



da; Im£o (ij ft, ix) Ym£' a (t, ft, ix) 

sinh 2 {ttx) Re 2 £rj (t, h, ix) + Im 2 £o (t, ft, ix) 



Evaluating the sum and int egra l close to T c q one can see 
that the first term in Eq. ( |C3|) is twice l arge r than the 
second one. Comparing Eq. (|C3|) to Eq. (C6) we obtain 

the old result: 3S<7 X 5 X +6 ^ = — 5<j x 3 x +4 \ ||, used later in 
Refs. |^, ||, and |l5| . But it is necessary to stress, that 
the last statement is not universal: far from the critical 
temperature, or at l ow te mper atures, close to H C 2 (0) , 
the integrals in Eqs. (C3) - (C6) are small with respect to 
the contributio n of t he su ms. Regarding the latter, they 
enter in Eqs. (C3) -(C6) with the opposite sign. After 
the summation in Sa x 5 x 8 ' these just cancel each other 

(in this region of temperatures 5a x 3 x +4 ^ ~ — Sai 5 ^ 6 ^). To 
avoid misunderstandingjid, it is more convenient to use 
the total contribution of the DOS-like diagrams 3-6 in 
the form: 
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4e 2 ft 



M 



dx Im 5m (t,h,ix) Im £ TO (i, ft., ix) 



7r 3 i 7_oo sinh 2 ttx Re 2 £ m (t, ft, ix) + Im 2 £ TO (t, ft,, ix) 



(C7) 



2. Asymptotic behavior 



c. Above the line H C 2 (T) but t <C ft C 2 (i) 



a. Vicinity ofT c0 , fields h « 1(# « # C 2 (0)) 



In this case lnt = e <C 1, the i/j-function in Eq.(Al4) 
can be expa nded. The function £ m (t, hAx) is determined 
by Eq.( |Al6| ). Its substitution to Eq.(C7) results in 



5a 



DOS 



14C(3). 



hi 



(l/2ft)-^(l/2 + ±) 



14C(3)e 2 J ln(l/e),ft<e 



(C8) 



This expression is valid in the vicinity of the critical 
perature T c q and exactly reproduces existing resultsl 



b. High temperatures, high fields 



Using the asymptotic Eq. (A22) one can perform the 
integration in Eq. (C7) and express 5a^® s in terms of 
the integral Jql : 



5a 



DOS 



7T 2 2^ 

m=0 



4/t(2ro+l) j n ih 



(2m + 1) 

Close to the line H C 2 (T) we can restrict ourselves to the 
LLL and immediately get 



Looking on the asymptotic behavior of J (r) at low 
temperatures, one notices that in contrast to the 
statement of Rcf Q the group of diagrams 5-8 does 
not give any contribution to 5a^£^ when temperature 
tends to zero. Nevertheless, a non-trivial contribu- 
tion of quantum fluctuations ~ In ft, found in Ref. [ [TBI , 
exists due to the regular MT term and diagrams 9 and 10. 



Next, we discuss the high-temperature asymptotic. As 
it was done above, we assume hat 1 and use Eqs. 
(gTj)-(|A2g). The sum over Landau levels in this case 
converges at large m max ~ t/h 1 and can be sub- 
stituted by an integral. The main integral contribution 
comes only from the region up to x ~ 1 and can be per- 
formed first. One gets 



5a 



DOS 

xx 



192 In t 



We see that this result differs from that on e of Rcf. 
The cancellation of the sums of Eqs. (|C3|) - (C6 in 5a, 



DOS 



removes the double logarithmic term lnlnt from it. Nev- 
ertheless, such terms in Saxx^ still appear from the reg- 
ular MT term and, as we will see below, from diagrams 
9 and 10. 

In the limit of high fields ft 3> t the summation over 
Landau levels gives: 



5a 



dos_ 7C(3U 2 e 2 ft\ 2 1 

384 IftJ m 22 4 



Appendix D: Renormalization of the diffusion 
coefficient: contribution of diagrams 7-10 

1. General expression 

We start with the calculation of diagram 7: 

QSK) = 2 e 2 T 2 ]T J -^L(q,Q k )X(q,e n ,Q k -e n ) 



(7) (7) 



(Dl) 



where the integrals of the Green's function products can 
be calculated in the standard way: 



^ (1) (^" n > £n-\-v 1 ^k—n) 



d D p 



>x (p) G (p, e n ) G (p, e n+u ) 



(2tt) d 

x G(q-p, n k - n ) = \kvvV<\ x t 2 6 (e n e n+v ) 9 (-e„O fc _ n ) 

(D2) 



r(7) _ r(7) (a x 

(2) — (1) \ k—n—v 1 iL k—ni c-n+v ) • 



Substitution of these expressions to Eq. ( pl| ) and ac- 
counting for the fact that £>q 2 = £>q 2 /2 results in 
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V xx -mv vei I ^ 2^ i, W> w fc ; ^ | 2e „ - n fc | + Vq 2 |2e„ + „ - fi fc | + Vq 2 \2e n + u „ - + Vq 2 



The 9— function 9 (e n £n+v) defines the limits of summation over fermionic frequencies n as (—00, —v — 1] and [0,oo). 
Changing the sign of summation in the first interval and then shifting the variable of summation e„ + lo v — > £„/ , one 
finds that the expression is even in Q^, which allows to present Q X J in the form of an analytical function of u> v , 
to perform the analytical continuation w v — > —iw, and to expend it over small u> : 

Q^ r (uj) = -8nP Ve 2 T 2 V [vq 2 L (q, Q k ) ^ V J , + — j) . (D3) 

W^ol^n + l^l+^q 2 ] 3 [2e n + \Q k \+Vq 2 ] 4 j 



The corresponding contribution to the c ond uctivity is de- 
termined by the imaginary part of Eq. ( |D3| ). Quantizing 
the motion of Cooper pairs and going over to the Landau 
representation, one finds 



k— — oc 



8£%(t,h,\k\) 
£ m (t,h, \k\) 
(D4) 

Comparing this formula with Eq. (B7) one can see that 
beyond the vicinity of T c q the contribution of diagrams 
7 and 8 given by the Eq.(D4) cancels the regular MT 
contribution [given by the first term of Eq. flBTp ]. 

Finally we proceed with the calculation of diagram 9. 
Two integrals of the three Green function blocks in it are 
equal and coincide with . Substituting Eq. ( [E>2| ) to 
the general expression for Qxxip^v) and performing the 
summation over fermionic frequencies in the spirit of the 
above calculations, one finds 



)(9) 



2M = - 



e 2 T 



00 p r i2 



[*i (|n fc | , Uv ) - * 2 (\n k+ »\,"»)} = Q§ + QgJ, 



where 
\l> 7 {x,ui v ) 



(D5) 



1 ^ + n x _ + w 

1 2 4ttT J Y \ 2 4ttT 



x + Vq 2 
4ttT 



(D6) 



with 7 = 1,2. 

There is no problem to per form analytical continuation 
of the first term of Eq. (D5): the function "fi (|fik| ,u) v ) 
is analytical in its argument oj v , and the corresponding 
contribution to Eq. ( p5| ) can be continued in the stan- 
dard way uj v — > —ioj — > 0. Expanding Eq. ( |D6| ) with 
7 = 1 over u> one finds the essential contribution to the 



electromagnetic response operator: 
Vv e 2 T 



Hi) 



3 (AttTY 
xL(q,n k )i>"' 



A;— — 00 



Vq 



2 _(Pq_ 

(27T) 2 

Vq 2 



4ttT 



(D7) 



The evaluation of the second term of Eq. ( D5 ) turns 
out to be much more sophisticated, since lo v appears in 
^2 (|^fc+j/| i w„) not only as parameter but also in the 
argument |f2fc_|_„| of this non-analytical function. The 
situation is analogous to the AL contribution and the 
same method of analytical continuation has to be ap- 
plied. The corresponding sum over bosonic frequencies 
is transformed in an integral over the contour C shown in 
Fig. [ll] with three regions of different analytic behavior: 



Q(2)K) 



2m 



Vq' 



j> coth^Z(g, 



d 2 q 
(2n) 2 

iz)ipi (|O fe+I/ | ,uj v ) . 



After shifting of the variable z of the integral over the 



line Imz = — lo v as — iz + w v — > —iz', one gets 
already as an analytical function of U3 V : 



)(9) 
1(2) 



Vv^e 2 f _ „ d 2 q 



Vq' 



dz coth 



(2n) 2 

x [V R (-iz + uj„,cj„)lmL R (q,-iz) 
+L A (q, -iz - uj v ) lva$> R (-iz, w v )] 



2TJ 



(D8) 



Obviously, this expression can be continued in uj v in the 
standard way uj u — > —iuj. 

We are interested in the imaginary part of Q^2 R {oS), 

i.e. only 1to^ r {— iz — icu, —iu) and Im^ r 2 i (— iz, —iu) are 
esse ntial. They can be can be written explicitly from Eq. 



lm^ R (—iz, -iu) = - 



3 (AttT)' 



-ReV>'' 



1 — iz 
2 + - 



Vq 2 



AttT 



(D9) 
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with Im*f ( 



ltd, 



5Im^f-(— iz, —itu). Since 



we are interested only in the linear w-part of ImQ^ ( 
in the analytically continued Eq. 



iui in the argument of L (q, 



(D8), one can omit 
Iuj) and recall that 



lmL A (q, 



-IvaL (q, —iz). One gets: 



which just gives a temperature independent constant. 
Let us stress that this constant is necessary for matching 
of the results in domain land Vllof Fig. fi. 



Vq 



4o>XV e 2 
3tt (AttT) 

1 —iz 

- A 

2 



x Re-0 



2 d 2 q 

(27T) 2 

h£>q 2 



4ttT 



dz coth ( 

\2TJ 



lmL R (q, —iz) . 



Now one can see that the integrand function is odd in 
z and its integration with symmetric limits gives zero. 
Hence, in linear approximation ImCjSuo;) = and the 
second term of Eq. (D5) docs not contribute to con- 



ductivity Going over to the dimensionless variables in 
Eq. (D7) and to the Landau representation, one finds 



that Sa^J = —Sa x 7 J/3. Finally, the total contribution of 
diagrams 7-10, determining the renormalization of the 
one-particle diffusion coefficient in the presence of super- 
conducting fluctuations, is 



8al n 



10 



4e^ 
3^6 



2 M 
m=0 fc= — oc 



(m 



1 8£%(t,h,\k\) 
2> £ m (t,h,\k\) ■ 

(D10) 



b. High temperatures, high fields 

In this domain of the phase diagr am, we cannot omit 
In t ^ 1 in the denominator of Eq. ( D1C ) , but the above 
consideration still is applicable. As a result we get 



Scr lx W (t > max{l, h}) = -^—r ( In In — In In t 



3tt 2 



In the limit of high fields ft> t 



T c0 t 



(D12) 



folx 10 ( h > max{l,t}) = ^1 ( In In J— - In In % 



37T 2 V T'cOT 



(D13) 



c. Above the line H C 2 (T) but t <^ h C 2 (t) 



2. Asymptotic behavior 



a. Vicinity ofT c0 , fields h < 1(H < H c2 (0)) 

In contrast to the AL, MT and DOS contributions, due 
t o p resence of the multiplier 2?q 2 in the num erator of Eq. 
( P5| ) [corresponding to (m+ |) in Eq. ( D10| ) close to the 

turns out to 



(DCR) 
(2) 



critical temperature T c o] , the value 8a 
be not singu lar in e at all. Substituting the summations 
in Eq. ( D1C ) by integrals, one finds 



8a 



.7-10 

xx 



( e « 1 ) = ^2 lnln 



" T c0 T 



0(e), (Dll) 



In this region one can restrict the consideration to the 
LLL approximation and use the asymptotic expression 
( |A22| ). In complete analogy with the case of the regular 
part of the MT contribution one finds 



Sal~ 10 {t<l,h) = ^lni. 



(D14) 



One can notice the close connection between the 8a 



and the 8a, 



MT(rcgl) 



7-10 

xx 



contributions, which is why the Eqs. 
( Dll| )-( D14 ) should be considered side by side with Eqs. 

(PFP¥ 
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